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In this talk I consider a two-dimensional lattice model for liquid crystals consisting 
of long rods interacting via purely hard core interactions, with two allowed orienta- 
^r) ■ tions defined by the underlying lattice. I report a rigorous proof of the existence of a 



nematic phase: by this I mean that at intermediate densities the system exhibits ori- 
entational order, either horizontal or vertical, but no positional order. The proof is 



based on a two-scales cluster expansion: first the system is coarse-grained on a scale 
comparable with the rods' length; then the resulting effective theory is re-expressed 
as a contours' model, which can be treated by Pirogov-Sinai methods. The talk is 
■ based on joint work with Margherita Disertori. 

CO 

o 

CO ; I. INTRODUCTION 

OV 

o: 

A nematic liquid crystal is a phase displayed by certain systems of anisotropic 
molecules, in which the distribution of orientations of the particles is anisotropic, 
^ 1 while the distribution of the particles in space is homogeneous. A picture to keep 

in mind is that of a system of cigar-shaped molecules in the three-dimensional 
continuum. The location of each particle in space is characterized by the position 
of the center of mass and by the direction of the main axis of symmetry of the 
molecule. In a nematic phase, all molecules' axes are directed approximately 
in the same way, while the centers are distributed in a random, "liquid-like" 
way. More formally, assuming the microscopic interaction among molecules to 
be translation and rotation invariant, in a nematic phase the underlying rotation 
symmetry of the system is completely broken, while the translation symmetry 
remains unbroken. 

Nematic liquids are physically characterized by peculiar refraction properties, 
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which can be conveniently used to filter polarized light, an effect that is exten- 
sively used in everyday's life in the liquid crystal displays of our computers and 
televisions. The technological importance of nematic liquids, which became clear 
already in the first half of the twentieth century, has been the boost for a large 
number of theoretical studies about the microscopic mechanism underlying their 
occurrence. Depending on the specific nature of the system under consideration, 
the intermolecular forces can be either of repulsive or attractive type. In certain 
suspensions of colloidal particles (e.g., aqueous suspensions of the tobacco mosaic 
virus) the dominant intermolecular interaction is a screened repulsive electrostatic 
force. 

For such systems, the theory of the microscopic mechanism underlying the 
spontaneous formation of a nematic phase is due to Lars Onsager 13]: roughly 
speaking, his idea was to picture each molecule as a long, thin rod, each char- 
acterized by a translational and a rotational degree of freedom. Since the inter- 
molecular interaction is repulsive and short ranged, it can be essentially thought 
of as a hard-core repulsion: therefore, the free energy is essentially equal to the 
entropy (i.e. the logarithm of the phase space volume corresponding to the al- 
lowed configurations of mutually non-overlapping molecules), which is the sum 
of a translational and a rotational entropy. At very low densities, so low that 
the average intermolecular distance is much larger than the rod's length £, every 
molecule has enough space to rotate freely around its center; therefore, the sys- 
tem is an isotropic gas phase. At higher densities, each molecule is surrounded by 
other molecules at an average distance smaller than I and, therefore, it is not free 
to rotate in all directions. In some cases, it may be favorable for the molecules 
to align spontaneously: it may be convenient for the system to substantially re- 
duce its orientational entropy, the loss being compensated by a much larger gain 
of translational entropy. To quantify this effect, Onsager computed an approx- 
imate equation of state by truncating at second order the virial expansion and 
by computing explicitly the first two Mayer's coefficients. Then he expressed the 
equation of state as a functional of the distribution function of the rods' orienta- 
tions and verified that for certain geometric shapes of the molecules there exists 
a non vanishing interval of densities where the free energy of the system is lowest 
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for an inhomogeneous distribution of the rods' orientations. 

Onsager's theory is the first example of an entropy- driven phase transition, 
i.e. an ordering transition induced by the competition of two entropic effects, 
each of which would separately like to make the system as disordered as possi- 
ble. The mechanism he proposed is often refereed to as the Onsager's excluded 
volume effect. Since then, there have been several attempts to prove in concrete 
models the correctness of this proposal, whose rigorous understanding is still an 
open problem in the context of three-dimensional systems with continuous trans- 
lational and rotational symmetries. On the other hand, there are indications of 
the validity of the Onsager's mechanism in the simpler context of long hard rods 
systems on Z d or on R d with discrete orientations. These can be thought of as toy 
models for lattice liquid crystals, in which the continuous rotational symmetry 
is replaced by a discrete one, much easier to treat; in this simplified context, we 
shall refer to a nematic-like phase (or nematic tout court) as a phase that breaks 
the discrete rotational symmetry, but not the translational one. In the following 
we shall briefly review, from a mathematical physics perspective, some important 
previous results concerning such systems. We shall restrict to systems where a 
hard core repulsion among rods is present and the set of allowed orientations is 
discrete; the case where the rod-rod interaction is of attractive nature and con- 
tinnons rotational symmetry is retained is another (very interesting) story, which 
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20] for some important results in this other 



we will not review here, see |2 
context. 



An outlook of previous results (mostly from a mathematical physics perspec- 
tive). 

1) In 1971, J. Lebowitz and G. Gallavotti [ijj proved the existence of orienta- 
tional order for a system of rods on a two-dimensional lattice. However, in their 
model orientational comes with positional order, which is not the case in nematic 
liquids. In this sense, the system studied by Lebowitz and Gallavotti is a model 
for a polarized crystal rather than for a liquid crystal. 



2) A few months after Lebowitz-Gallavotti's paper, D. Ruelle 19( succeeded in ex 



tending the use of the Peierls' method to a system of infinitely thin rods in the 2D 
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continuum (Widom-Rowlinson model). Such model exhibits orientational with- 
out positional order. In this respect, it represents the first rigorous confirmation 
of the Onsager's excluded volume effect in a system with discrete orientations. 
However, the assumption of infinitesimal size is unphysical and drastically simpli- 
fies the analysis: horizontal rods interact only with horizontal rods, and similarly 
for vertical. It is then natural to ask whether the same result remains valid for 
finite-size molecules. 

3) The simplest model to look at for studying finite size molecules is the k-mer 
model on Z 2 : particles are modelled as 1 x k rectangles placed on the sites of Z 2 , 
either horizontally or vertically, interacting via a hard core repulsion; the density 
of /c-mers is fixed and smaller than the close packing density (i.e., the rods do 
not cover the plane, but leave a prefixed fraction of the total volume empty). I 



k = 2 we have dimers at finite density. In 1972, O. Heilmann and E. Lieb 11] 
proved that the thermodynamic and correlation functions of dimers on Z 2 are 
analytic for all finite densities. Therefore, hard-core dimers do not exhibit any 
nematic-like phase. 

4) In an attempt of constructing a system of finite size anisotropic molecules ex- 
hibiting a nematic-like behavior, in 1978, O. Heilmann and E. Lieb 12| proposed 
a variant of the dimer system, by adding attractive forces between aligned dimers. 
By reflection positivity, they proved the existence of orientational order. Presum- 
ably, orientational order comes without translational one 16], as it should in a 
nematic phase. However, the attractive forces completely change the mechanism 
driving the ordering transition. We are still left with the question of whether 
the excluded volume effect enough to induce order. What about fc-mers, with 
k > 3? To put in in Heilmann and Lieb's words: It is doubtful whether hard rods 
on a cubic lattice without any additional interaction do indeed undergo a phase 
transition. 



5) In 2006, D. Ioffe, Y. Velenik and M. Zahradnik [13] proposed a "polydisperse" 
variant of the k-mer model, where rods of all possible lengths are allowed, with 
statistical weight depending on k. They mapped the model into the 2D Ising 
model and by the exact solution they proved the existence of an isotropic-nematic 
transition, which is the image under the mapping of the usual paramagnetic- 
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ferromagnetic transition in the Ising model. This result strongly suggests that 
the same transition should take place in the pure k-mer system, but unfortunately 
the method of [3] breaks down in the presence of apparently harmless changes 
in the relative rods' weights. 

6) In 2007, A. Ghosh and D. Dhar |9] numerically identified in a very clear man- 
ner a nematic phase for pure k-mers with k > 7 at intermediate densities. Their 
result further confirmed the expectation arising from the work of [ljj. 



II. THE MODEL AND THE MAIN RESULTS 

Motivated by these results, M. Disertori and I recently reconsidered the k-mer 

n 

problem on the two-dimensional square lattice [5J . An informal statement of our 
main theorem is the following. Let p G (0, 1/k) be the density of rods. 

For k large enough, if k~ 2 <C p <C k~ l , the system admits two distinct infinite 
volume Gibbs states, characterized by long range orientational order (either hori- 
zontal or vertical) and no translational order, selected by the boundary conditions. 

This result is the first rigorous proof of the existence of a nematic phase (i.e., a 
phase that breaks rotational symmetry but does not break translational symme- 
try) and of the validity of the Onsager's volume effect in a system of anisotropic 
particles of finite fixed size. It solves a conjecture posed more than three decades 
ago, as explained in the historical outline above. Before I give a sketch of the 
proof, let me give a more formal statement of the same result, which makes clear 
how the boundary conditions are chosen, and clarifies certain key properties of 
the limiting Gibbs states. Even if more formal, the definitions and statements 
below will still be slightly unprecise, in order to avoid cumbersome notations or 
digressions; for a mathematically precise exposition, I refer to {5]. 

Consider a box A C Z 2 of side L (eventually, L — > 00). We denote by the 
set of allowed rod configurations in A (we shall say that a rod belongs to A if its 
geometrical center do so; similarly, a rods' configuration belongs to A if all its 
rods belong to A). Let also (resp. Q^) be the set of horizontal (resp. vertical) 
allowed rod configurations in A. We define the grand canonical partition function 
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with open boundary conditions as: 

Z (A) =J2 zlRl 
Ren A 

where z is the rods' activity (independent of the orientation) and \R\ is the number 
of rods in the configuration R. Similarly, the partition function with + boundary 
conditions is: 

Z(A\+) = zW 

A 

where P™ 1 is the internal peel of A, i.e., a 2/c-thick boundary layer of A. In 
other words, the partition sum with + boundary conditions is such that all the 
rods with centers closer than 2k to the boundary are forced to be horizontal, 
while those in the "bulk" have no constraint on their orientation. An analogous 
definition is valid for — boundary conditions, with the role of the horizontal rods 
exchanged with the vertical. 

In addition to the partition functions, we can define as usual the correlation 
functions with open, or +, or — boundary conditions. For example, the correlation 
functions with + boundary conditions are defined as follows. Let Ax be a local 
observable with finite support X, then: 

RG^Ia '■ 

Rlpint G^p int 

and (Ax) + = lim (Ax)t- The collection of the limiting statistical averages 

A|— >oo 

(Ax)X as Ax spans a complete set of observables defines the infinite volume 
Gibbs state with + boundary conditions, (-)X- Analogous definitions hold for 
(■}~ or (-) . In terms of these definitions, a more precise way of formulating our 
main result is the following. 

r-i 

Theorem [5|]. There exist Kq, Eq > suchthat, ifk > Kq andm&x{zk,e~ zk2 } < 
Eo, then the two infinite volume states (•) ± exist. They are translationally invari- 
ant, puref? } and are different among each other. In particular: 
1) If E7 is the event that the rods in a window A^ of side I = k/2 centered at 
£o are all vertical, then 

Prob+(E-) := (x(E7)) + < e~^ , 
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where x{E^ ) is the characteristic function of the event and c is a suitable positive 
constant. 

2) If n x is the indicator function that is equal to 1 if a rod has a center in x G Z 2 
and otherwise, then 

p=(n x ) + = (n x )~ = z(l + 0(e c ')) 

pip -y) = (n x n y ) + = (n x n y )~ = p 2 (l + 0(e c '^ k ^ 
where e = maxjzfc, e~ zk2 } and d is a suitable positive constant. 

III. IDEAS OF THE PROOF 

While we refer to jlj] for a complete proof, we sketch here the main ideas 
involved in the proof, and its main steps, just to give a flavor of the methods 
involved. 

1) The first step consists in coarse graining A in square tiles A of side t = k/2: 
this means that we partition the original box in a collection of tiles, each of 
which contains in average many (~ zk 2 ^> 1) rods; in this respect the tiles can 
be thought of as being mesoscopic. On the other hand, each tile A is so small 
(the side is half the rods' length) that only rods of the same orientation can have 
centers in A, due to the hard rod condition. 

2) Given a tile A, once we prescribe the orientation of the rods with centers in it, 
the effective interaction between rods of the same orientation, say horizontal, is 
weak: in fact, the hard core repulsion just prevents two rods to occupy the same 
row, an event that is very rare, since the density of occupied rows is ~ 4p <C 1. 
It is remarkable that the standard cluster expansion allows us to quantify how 
close to Poissonian is the distribution of centers in A, once the orientation is 
prescribed. 

3) Each tile can thus be of three types: +1 (horizontal), —1 (vertical) and 
(empty). In this way, we can associate with every allowed rods' configuration 
on A a corresponding spin configuration on the coarser lattice of the tiles' cen- 
ters. By summing over all the rods' configuration corresponding to a given spin 
configuration, we are left with a partition sum over spins, which defines an effec- 
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tive spin model. This spin system has the following features: (i) the interaction 
between spins is short range and strongly attractive, due to the hard core; (ii) 
The vacuum configurations are unlikely, since the probability that A is empty is 
^ e -(const.)zfc 2 Therefore, the typical spin configurations consist of big connected 
clusters of "uniformly magnetized spins" , separated by boundary layers (the con- 
tours), which contain zeros or pairs of neighboring opposite spins. 
4) We are left with studying this effective contour theory. The idea is to use a 
Peierls argument or cluster expansion methods. However, there are a few issues 
that make life complicated. 

(4. a) First of all, the inter-contour interaction turns out to be an exponentially 
decaying iV-body interaction, with N arbitrary; if the nature of this iV-body in- 
teraction were generic then we would really be in trouble: nobody knows how to 
work out a convergent cluster expansion for generic iV-body interactions, even 
if exponentially decaying. Luckily enough, the iV-body interaction we need to 
deal with is quasi- ID, in the sense that only contours whose horizontal or vertical 
projections overlap interact among each other. This makes the problem treat- 
able. The right strategy is to follow Brydges' suggestion |3|] : if at first you do 
not succeed, then expand and expand again! In other words, we perform a sec- 
ond Mayer expansion of the multi-contour interaction and we collect together the 
connected components (polymers); the resulting polymers have purely hard core 
interactions, and the polymers' activities satisfy similar bounds as those of the 
original contours. 

(4.b) The activity of the contour is defined in terms of a ratio of partition func- 
tions that at the beginning of the story we do not know how to compute (yet). 
Moreover, the contour theory is not exactly symmetric under spin flip, due to 
the finite size of the rods (i.e., the activity of a given contour is not the same 
in the presence of + or — boundary conditions, unless the shape of the contour 
itself is rotation invariant). To solve these two issues, we use the methods of 



Pirogov- Sinai theory 



18 



21] . By induction, we show that the polymers satisfy a 
Peierls' condition, i.e., the probability that a given contour or polymer occurs is 
exponentially small in the size of its geometric support. For details, see jlj]. 
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IV. CONCLUSIONS AND OPEN PROBLEMS 



I reported the first proof of the existence of a nematic-like phase in the pure k- 
mer system in two dimensions, for k sufficiently large and intermediate densities. 
This proof solves a long standing conjecture about the validity of the Onsager's 
excluded volume effect in systems of finite size anisotropic molecules. Of course, 
many questions remain open. Let us mention a few of them. 

1) Three dimensions. It would be nice to extend the results of jlj] to three dimen- 
sions. Presumably, an analogous result should remain valid for 1 x 1 x k rods 
on Z 3 . A more interesting question is what happens for anisotropic molecules of 
the form 1 x fc* x fc (0 < a < 1) with six allowed orientations and centers in M 3 : 
in this case a priori there is the possibility that other kinds of liquid crystalline 
order (smectic, chiral, etc 4|) emerge. 

2) High density. As p is increased towards close-packing, a second transition to 
a "disordered" state (i.e., a state with no orientational order) is expected {9}. 
The understanding of this dense phase is completely open. In particular, at close 
packing, the best known variational bound on the entropy is: > C s ^§^ for 
an explicitly known constant C s , which is obtained by making a "striped" ansatz 
Q]. A very interesting open problem is to prove that & ~ asymptotically 
as k — > oo. 

3) Rotational invariance. "Of course", the most important and difficult open 
problem is to prove the existence of orientational order in a model with gen- 
uine rotational invariance. There are some examples of two or three-dimensional 
gases with a continuous internal degree of freedom, interacting via a tensorial 
attractive interaction, where the existence of nematic order (or quasi-order in 



2D) can 



3e proved by a combination of cluster expansion and reflection positivity 



20| . However, the understanding of the Onsager's excluded volume effect 



is completely open in this case. I think that this last problem is substantially 
more difficult than the previous two. We have a serious lack of understanding of 
continuous symmetry breaking phenomena, as witnessed by the state of art in the 
problem of existence of spontaneous magnetization in 3D classical and quantum 
ferromagnets: there are only very few special cases where spontaneous magnetiza- 
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tion can be rigorously proved , the proof being based on reflection positivity, 
which is very fragile under apparently harmless changes in the Hamiltonian. It 
is very likely that future progress on the theory of the Heisenberg ferromagnet 
will also help us in understanding the validity of the Onsager's excluded volume 
effect in continuous three dimensional systems. It may also be of help to try to 
attack the issue of continuous rotational symmetry by first looking at mean-field 



like models, possibly in the spirit of 
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